Fluid Velocity Fluctuations in a Suspension of Swimming Protists 
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In dilute suspensions of swimming microorganisms the local fluid velocity is a random superposi- 
tion of the flow fields set up by the individual organisms, which in turn have multipole contributions 
decaying as inverse powers of distance from the organism. Here we show that the conditions under 
which the central limit theorem guarantees a Gaussian probability distribution function of velocities 
are satisfied when the leading force singularity is a Stokeslet, but are not when it is any higher 
multipole. These results are confirmed by numerical studies and by experiments on suspensions of 
the alga Volvox carteri, which show that deviations from Gaussianity arise from near-field effects. 
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A key feature of the inertialess world inhabited by 
microscopic organisms is the very long-range flow fields 
they create as they swim. For neutrally-buoyant, self- 
propelled organisms the far-field behavior of the velocity 
is that of the force dipole (stresslet) created by the op- 
posed actions of their flagella and cell body on the fluid. 
Theories incorporating such fields in the fluid stress ten- 
sor [1L and simulations of suspensions of dipolar organ- 
isms [2J have shown the formation of large coherent struc- 
tures that are highly suggestive of those seen in exper- 
iments on the bacterium B. subtilis Q. The suggestion 
Q that hydrodynamic interactions underlie these vor- 
tices and jets was made by analogy with the appearance 
of similar patterns in suspensions of sedimenting parti- 
cles 0| , although interactions in the latter are due to the 
force monopole (Stokeslet) fields arising from the den- 
sity mismatch between the particles and fluid. Although 
Stokeslet and stresslet fields have different topologies, it 
is striking that the two systems display similar coherence. 

The relationship between suspensions of microorgan- 
isms and of sedimenting particles takes on new signifi- 
cance in light of measurements of velocity fields around 
freely-swimming organisms which emphasized that 
the Stokeslet field dominates that of the stresslet beyond 
a length A ~ Td/F g , where d is the offset between the 
flagellar thrust T and the body drag, and F g is the net 
gravitational force per organism. A can be surprisingly 
small when compared to the organism radius R: while 
for the unicellular alga Chlamydomonas reinhardtii [6[ 
(R ~ 5 //m) A ~ 30R, for its multicellular descendant 
Volvox carteri [7] (R ~ 200 — 400 /xm) there is the strik- 
ing conclusion that A ~ R; the Stokeslet dominates the 
flow field. It was therefore suggested Q that suspensions 
of Volvox would be more similar to those of sedimenting 
particles than previously thought, the chief difference be- 
ing the component of the organism's motion from active 
swimming. Hence there is fundamental interest in the 
question: What are the statistics of fluid velocity fluctu- 
ations in a suspension of swimming microorganisms ? 



Here we present theory, experiments and simulations 
that elucidate a number of aspects of this question. We 
determine the condition on the leading force singularity 
of a swimmer in order that a random superposition of 
its velocity field has a Gaussian probability distribution 
function (PDF). The condition admits the Stokeslet but 
excludes the stresslet and higher multipoles, so the pres- 
ence or absence of density matching has a qualitative 
effect on the statistics. The velocity distribution func- 
tions found in experiment and simulation display clear 
non-Gaussian tails which we suggest arise from near-field 
effects dJ. The large size of Volvox allows study of the 
scale of fluctuations as a function of the number of organ- 
isms at fixed container size, complementary to the lim- 
iting procedure often adopted in sedimentation [§]. Our 
result complement recent studies of the short-time PDFs 
of tracer particles in suspensions of Chlamydomonas [Io| , 
where non-Gaussianity was fou nd, and to studies of fluc- 
tuations in bacterial baths 



11 



Consider a suspension in a box of linear scale L, with 
N swimmers of radius R. If the volume fraction (f> — 
4irR 3 N/3L 3 is sufficiently small, the PDF of velocities 
reflects the statistics of a random superposition of the 
flow fields around each swimmer. For a uniform spatial 
distribution of swimmers, averaging over their positions 
is equivalent to integrating over space with the swimmer 
at the origin. If the velocity around a swimmer decays as 
|v(r)| ~ A(fl)/r n , with ft standing for angular variables, 
the probability distribution P(v) of velocities is 
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assuming a spherical container. The tail of the distri- 
bution can be determined from the behavior of P un- 
der the rescaling r — > ar. Since S(v — A/(ar) n ) — 
a n 5(va n — A/r n ), and noting that for large v the ar- 
gument of the 5-function vanishes at small r, we deduce 
that the integral does not depend on the upper limit L/a 
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FIG. 1: (color online). Experimental setup and measured velocity fluctuations, (a) Schematic of the imaging and illumination 
system, (b) Experimental ^-component (and x-component, inset) of the fluid velocity in a Volvox suspension (N = 300 cm -3 ) 
as a function of time in the central PIV grid domain of the chamber. Red lines indicate velocities in SVM, while blue lines are 
with approximate density matching of the external fluid: SVM+3% v/v Percoll (Sigma). 



(which can be taken to oo), and hence 
P(v) = a 3+n P(va n ) -> P(v) oc 
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The second moment of P{v) is finite only if n < 3/2, the 
case of a Stokeslet (n = 1). This is the condition for 
validity of the central limit theorem; the velocity field 
from a large number of independently placed Stokeslets 
is Gaussian. It will not be so for any higher integer sin- 
gularity, such as stresslets (n = 2) or source doublets 
(here termed 'sourcelets') (n = 3) 13]. If below a certain 
radius the decay law deviates from v oc r~ n , the PDF 
shape ([2]) will break down above the corresponding value 
of v. Hence, deviations from Gaussianity provide a direct 
probe of the near-field velocity around the swimmers. 

The spherical colonial alga Volvox is a remarkably use- 
ful system for the study of many aspects of biological fluid 
dynamics 13- 17] because of its size, high symmetry, ease 
of growth, well-characterized biology, and the existence 
of a range of mutants. In our experiments, Volvox car- 
ters f. nagariensis strain EVE were grown axenically in 
SVM [lH in a diurnal growth chamber set to a cycle of 16 
hours artificial cool daylight (~ 4000 lux) at 28° C and 8 h 
in the dark at 26° C. We used synchronized colonies from 
the first day of the 48 hour life cycle to obtain the highest 
motility. A concentration c = 10 — 500 cm~ 3 (a volume 
fraction below = 0.015) of organisms was prepared in 
SVM, with added 2 fj,m polystyrene seeding particles or 
6 jUin tracer particles (Polysciences) at a concentration 
of ~ 25 ppm, and placed into a glass container (lxlxl 
cm). The container was plac ed a thermal bath (Fig. QJi) 
to eliminate convection [19|, and was illuminated with 
a thin laser sheet (< 300 /im) from a 100 mW, 635 nm 
laser (BWTEK). Video was captured at frame rates of 
0.4 - 5 Hz by a CCD camera (Pike F145B, Allied Vi- 
sion Technologies) connected to a long-working distance 



microscope (Infinivar CFM-2/S, Infinity Photo-Optical). 
The fluid velocity was measured by PIV (Dantec Dynam- 
ics), typically producing a 63 x 63 rectangular lattice of 
velocity vectors. Alternatively, we measured tracer and 
Volvox trajectories by custom (Matlab) PTV software. 

Our simulations of protist suspensions used a model in 
which the velocity field created by a Volvox is the sum of 
a downward-pointing Stokeslet and a sourcelet, 
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where n is a unit vector along the colonial axis. The 
sourcelet represents the near-field flow found by direct 
measurements [B| and in a model with a constant force 
density distributed over the colony surface [l5[ , and is im- 
portant for the statistics of high fluid velocities. Both sin- 
gularities are cut off at the colony radius. The Stokeslet 
strength was fixed by an empirical fit to data on the sed- 
imentation velocity as a function of R [l6| [v sc d — uR, 
with a = 1 s _1 ], while the relative sourcelet strength 
ji was studied in the range < fi < 10. We consider 
the motion of colonies within a non-interacting 'ideal gas' 
model [20] which, despite its simplicity, gives satisfactory 
agreement with the experiment; the fluid is unbounded, 
the swimmers are confined to a rectangular container 
(cage) with reflecting walls, and the position Xj of the 
jth swimmer and its axis vector evolve as 
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where v p is the propulsion velocity, Wj and Wj are white 
noises with diffusion constants D and D, in 3D and on 
a unit sphere, respectively. They represent the random 
influences on the motion of Volvox - irregularities of flag- 
ellar beating and, partially, mutual advection of colonies. 




FIG. 2: (color online). Statistics of velocity fluctuations. Colored circles are experimental data for a suspension of Volvox 
with mean diameter 300 /jm. Colored triangles are corresponding numerical simulations excluding sourcelets; solid lines at 
right are for simulations including sourcelets (fi = 4). Individual colors indicate different numbers of colonies in container: 
green (11), red (42), and blue (210). (a) PDF of fluctuations in horizontal velocity v x . (b) Standard deviation of fluid velocity 
normalized by colony sedimentation speed u se d [z-component (squares), s-component (triangles)], for colony mean diameter 
460 jUm (open symbols), and 220 /im (solid symbols), (c) Central region of PDFs of v x normalized by their standard deviations, 
for N = 11, 28,42, 128 and 210. Dashed black line is a Gaussian fit. Full PDFs of (d) v x and (e) vertical velocity v z . 



The latter is negligible for the most part, since the typ- 
ical velocity of the resulting flow is found to be much 
smaller than v p . This is not true, however, when two 
or more colonies come close. Although such events are 
relatively rare, they are important for uniformizing the 
spatial distribution of Volvox: without them the bottom- 
heavy colonies would gather in the upper part of the con- 
tainer, contrary to observations. For the same reason, in- 
cluding the bottom-heaviness and the sedimentation into 
([3]) in the absence of mutual advection would be inconsis- 
tent. The primary (and minor) consequence of neglecting 
bottom-heaviness is this model does not reproduce the 
angular distribution of the colonies' axes. Inclusion of 
sedimentation makes only minor changes to the results. 

An example of experimental measurements is the time 
trace of local fluid velocity in the center of the sample 
chamber (Fig. [Tp). We see that the observed fluid motion 
is created primarily by the Stokeslets of the swimmers, 
for when the fluid density was increased to match the 
density of Volvox, the typical velocity fluctuations were 
reduced drastically. Yet, the peaks due to the near-field 
source doublets (from a swimmer passing very close to 
the observation point) remained undiminished. 

On a more quantitative level, we examined the PDF of 
velocity fluctuations (Fig. [2J as a function of the number 
of colonies in the container, and at various stages in the 
lifecycle, so that the colony size and sedimentation speed 
vary over a significant range. Data for the smallest num- 
ber of swimmers in the container shows a clear power- 
law tail consistent with the form v~ 4 expected from Eq. 
[2j and in agreement with simulations done with pure 
Stokeslets. As expected from a gas of Stokeslets, the PDF 
of the velocity shows convergence to a Gaussian with the 



number of swimmers: for 210 swimmers the Gaussian- 
ity persists up to 2.5 standard deviations (Fig. [2t), but 
with clear tails (discussed below). Once normalized by 
the sedimentation speed, the standard deviation of the 
velocity collapses, showing that the fluctuations are pro- 
portional to the Stokeslet strength (Fig. Wp)- In an ideal 
gas of Stokeslets, the standard deviation of the veloc- 
ity fluctuations grows as cx >/~N by virtue of the central 
limit theorem. In the presence of swimmer correlations 
it should grow faster, but no faster than cx N. The ob- 
served law lies between these two powers, much closer to 
the former (Fig. O)), supporting the ideal gas approxi- 
mation, and distinct from the result TV 1 / 3 found in sed- 
imentation [4J, where the mutual advection of particles 
in each other's Stokeslet fields is the sole contribution to 
velocity fluctuations. Fluctuations in Volvox suspensions 
are stronger for larger swimmers, due to their larger sed- 
imentation velocity (stronger Stokeslets), and the ratio 
1Z = o-(v z )/a(v x ) is close to 2 for all N. This is found in 
the numerics with Stokeslets+sourcelets (whose orienta- 
tions are uniformly distributed). Without sourcelets, the 
numerics yield R~2.8k 2y/2. For a single Stokeslet off- 
set by a sourcelet, the ratio can be computed analytically, 
averaging over the swimmer's position being replaced by 
spatial integration, as in ([1]); 1Z ranges from 1 (for a ran- 
domly directed sourcelet) to 2^2 for a pure stokeslet. 

Inclusion of a sourcelet in the simulations results in 
tails in the PDFs similar to the experimental data (Fig. 2 
right). This allows us to conclude that the observed tails 
are due to the near-field component of the swimmers' 
flow. The tails appear to be exponential, but the range 
of our data is insufficient to prove this. For example, 
the tail of the data for N — 210 is equally well fit by 
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The greatest challenge is a theoretical understanding of 
the form of the non-Gaussianity, which is known to ap- 
pear as well in other contexts, such as inelastic gases [2l|. 
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FIG. 3: (color online). Power-spectra of velocities. Solid 
lines and dotted lines are experimental spectra of v z and v x , 
respectively, for suspensions with mean Volvox diameter 290 
/im, for N = 20 (red), 70 (green), 250 (cyan), and 406 (ma- 
genta). Solid and dotted black lines are numerical results 
for N = 210 (rescaled in y), with ^ = 4. Inset: collapse of 
power-spectra in numerics with varying container sizes (0.5-2 
cm) and v p (150-600 fj,m/s), with it, from 16.7s to 66.7 s. 

P oc (v 3 + Nv 3 )~ 1 with = 25 /im/sec. A similar 
situation occurs in the PDF of the vertical velocity, where 
the core convergence to a Gaussian is less advanced due 
to the inherent asymmetry of the Stokeslet field. 

The velocity power-spectra show a decay close to f~ 2 
(Fig. [3]), suggesting a Lorentzian power-spectrum of fluc- 
tuations P(f) = (f 2 + (27tt)~ 2 ) -1 , i.e. an exponential 
velocity autocorrelation (v(O)v(t)) cx exp(— t/r). Sup- 
plemented by the Gaussian PDF, this equation amounts 
to modeling the velocity fluctuations as an Orstein- 
Uhlenbeck stochastic process. The motion of Volvox is 
primarily deterministic. For a concentration c ~ 100 
cm" 3 the mean free path can be estimated as (ttR 2 c)~ 1 ~ 
10 cm, which is larger than the container size of L = 1 
cm. Thus, the deterministic term in (|3]) sets a ballistic 
time Tb = L/vp ~ 30 s which is smaller than the diffu- 
sive time scale = L 2 / D ~ 100 s or the dephasing time 
Tph = 1/-D ~ 100 s. Hence, it is t\> that sets the corre- 
lation time in this ideal gas model. We checked it in the 
numerics (Fig. Q^inset)), and the characteristic r in the 
experimental data is close to that in the numerics. 

In summary, we have introduced a connection be- 
tween the statistics of velocity fluctuations in suspen- 
sions of swimming protists and the type of force singu- 
larity associated with the organism motion. Experiments 
and numerical results show clearly the existence of non- 
Gaussianity in the velocity PDFs, which are suggested 
to arise from the details of fluid flow near the organisms. 
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